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ABSTRACT 


Land use land cover change (LULCC) processes are directly or indirectly product of underlying 
causative factors playing out on the landscape. This study examines the contributing factors 
responsible for LULCC in Jos Plateau state using two different methods. The Analytical 
hierarchical process (AHP) and Binary models were used to prepare sixteen (16) suspected 
external driving factors (EDF's). Euclidean distance analysis was carried out on the proximity 
EDFs as well as reclassification based on AHP scale I — 9 and binary mode as I & 0. Saaty’s 
pairwise comparison matrix was employed to generate the weights for all the EDFs with a 
consistency index (CI) of 0.1 7was achieved. Suitability images were generated by multiplying each 
EDF with the corresponding weight. The fuzzy set membership standardization was carried out on 
the AHP image through sigmoidal function type with monotonically increasing membership 
function.The multiple regression technique was used to measure the sensitivity of each driver 
against the change period (1986 — 2019). The result shows that five (5) EDFs where consistent in 
both approaches (AHP & Binary). They are: distance to conflict areas, distance to major roads, 
LULC, distance to settlements and distance to river/ channels. It is therefore, appropriate to 
examine EDFs to know how each driver is contributing to the changing landscape in the study 
area. The findings has revealed the main factors driving LULC changes in the study area and has 
provided a reference frame for assisting in the development of sustainable land management and 
ecological protection policy making and decisions. 
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1.0. Introduction 


Sensitivity investigation of external driving factors (EDFs) is a pointer that estimates the variation in 
the model output with respect to change in the model’s input parameters. It is also used for estimating 
the strength of the model. Burg (2016) asserts that feeding a large amount of data into the model does 
not always result in better outputs. In fact, the more the data fed into the model, the higher the chances 
of introducing uncertainty and hence the more error. In this research, sensitivity analysis was 
performed because the drivers are products of varying processes and prepared at different resolutions. 
Many of these terms exists, most commonly used are: anthropogenic driving factors in Yichen et al. 
(2014), influencing factors by Maher et al. (2016) & (2017). Marcelo et al. (2008) referred to it as 
proximate causes/ underlying driving factors and Kim et al. (2014) maintained as driving factors; 
Jiujun (2017) as human factors while in Sanchayeeta et al. (2017) sees it as proximate/ indirect causes 
of LULCC. Land use land cover changes (LULCC) are significant changes driven by human activities 
(Alemayehu, 2019). In the history of man, the land has been tightly attached to economic, social, 
infrastructure and other human activities (Lambin et al., 2003). These changes affects life support 
functions and humans livelihood (Lambin et al, 2001; Lambin and Geist, 2006). Significant LULCC 
is caused by a combination of varying factors depending on each locality condition (Hassen et al., 
2015; Eyayu et al., 2009; Woldeamlak and Sterk, 2005; Eleni et al., 2013; Mohammed, 2011; 
Woldeamlak, 2002; Hazem, 2020). 
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As observed in the literatures above, different terms are used but they all refer to the same thing. Even 
though, each term may be given different interpretation under different scenario. This has set in some 
lack of clarity when dealing with some of these driving factors of LULCC. Therefore, this research 
adopted the term ‘External Driving Factor (EDF)’for consistency and clarity. Various EDFs as used 
by different studies and understanding the underlying factor (driving mechanisms) of LULCC caused 
by a variety of driving forces is one of the optimal goals of global change research in recent decades 
(Xiangmei et al., 2016). To understand the human and biophysical processes of LULCC, many 
researchers focus on the various forces driving LULCC, including socioeconomic by Xie et al. 
(2005), demographic in Shi et al. (2010), political by Kanianska et al. (2014), technological as in 
Hasselmann et al. (2010), biophysical by Serra et al. (2008) and industrial structure in Shu et al. 
(2014) to provide effective support for developing urban land planning and management regulations. 
To comprehensively analyse the driving factor’s effects and mitigate the negative impacts of LULCC, 
Shu et al. (2014) investigated the effects of various factors, including natural Eco environment 
factors, land control policies, accessibility factors, and neighbourhood factors on urban land 
expansion during various periods in different regions. Chen et al. (2014) selected industrial structure, 
Gross Domestic Product (GDP), transportation, and policy as the driving factors to study the impacts 
on urban land expansion and sustainable urban development. Bandit (2018) analysed the spatial 
relation between independent and dependent variables using Multi binary logistic regression 
technique. The study concluded that, both biophysical and socio-economic variables had meaningfully 
contributed to LULC conversions. 


Attempts have been on analysing these contributing drivers to LULC dynamics in several parts of the 
world, notable are Pena et al. (2007); Ellis (2010); Alemayehu et al. (2019) who carried out a 
comprehensive studies of landscape change to understand the underlying processes and a full range of 
approaches from the natural and social sciences. Logistic regression model was used to establish the 
association between socio-economic drivers and LCC in (Dimobe et al., 2015). Similarly, Yuliana and 
Kaswanto (2016) applied the linear regression equation using the stepwise regression to analyzed 
drivers of LULCC. Sanchayeete et al. (2017) identified influencing factors of deforestation and 
simultaneous plantation driven reforestation in Bannerghatta National Park using binary logistic 
regression likewise the same technique was applied in Mondal et al. (2014) to investigate and 
analyzed the drivers of LULCC. Due to the fact that a correlation may exist among the selected 
variables, principal component analysis (PCA) and general linear model (GLM) were employed to 
identify the relative significance of the driving factors (Xindong et al., 2014). Geospatial analysis 
method was used to identify drivers of LULCC in (Janina et al., 2017). Yichen et al. (2014) carried 
out a qualitative analysis of the impacts of drivers on LULC change using Pearson’s correlation 
coefficient where correlation among the built-up areas and the socioeconomic statistics was analyzed. 
Briassoulis (2000, 2008) emphasized the reason why the linkage between LU and LC change in 
environmental impacts and their contribution to global change are mediated to a considerable extent. 
Thus, their analysis needs the examination of the ways in which LU relates to LC change at various 
levels of spatial and temporal detail. LULCC is usually categorized into two broad classes: conversion 
(a change from LULC category to another e.g. from forest to grassland) and modification (a change 
within one LULC category e.g. from rain fed cultivated area to irrigated cultivated area) (European 
Commission, 2001). LC conversion is the complete replacement of one cover type by another such as 
deforestation to create cropland or pasture (Rachmad, 2014) and these are responsible by some 
underlying factors. Hassen and Muhamed (2017) identified some major causes of LULCC in Tana 
watershed. 


Uncontrolled landscape expansion is a fundamental problem experienced in many societies of the 
world today as a result of population explosion. This has been the leading causes of LULCC today in 
so many parts of the world. Jos metropolis is expanding in all directions (Adzandeh et al., 2015) 
creating serious environmental problems. Therefore, this research seeks to analyse these LULCC 
driving factors (external driving factors) to measure how sensitive they are in contributing to the 
present landscape change in the study area. 
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2.0. Methodology 


2.1. Study area 


Jos-Bukuru metropolis the study area lies within Jos North and Jos South LGAs of Plateau state. 
Plateau state as a whole derives her name from Jos plateau (geographical landscape) that dominates 
this part of country and is located relatively at the Centre of the country Nigeria. Jos lies between 
longitude 8° 40' 47.15"E and 9° 8' 13.20"E, latitude 9° 36' 11.54"N and 10° 2' 20.51"N (Figure 1) 
covering an area of about 7,780 sqkm with a population close to 900,000 and a population density of 
391 persons per km* making it to be the most densely populated part of plateau state. 


MAP OF NIGERIA SHOWING PLATEAU STATE 
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Figure 1: Map of study area 


2.2. Materials 


Table 1 shows the materials used for this research. 


Table 1: Materials used 
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2.3. Methods 


2.3.1. Determination of external driving factors (EDFS) of LULCC 

Deriving external factors to be used for further analysis began by first identifying suspected factors as 
put forward by some authors Al-sharif and Pradhan (2014); Shafizadeh-Moghadam and Helbich 
(2015); Arsanjani et al. (2013); Chowdhury and Maithani (2014); Zeng et al. (2015) and Maher et al. 
(2017). Some of the factors put forward by these authors were adopted and combined with local 
factors peculiar to the study area as distance to conflict areas and distance to mining areas. In this 
research, a total of sixteen (16) suspected EDFs (variables) having spatial characteristics were 
identified and grouped into two as follows: (1) geophysical and (11) proximity as shown in Table 1. 
Spatial thresholds were adopted for these layers created as shown in Table 2. They were then 
weighted and ranked using Saaty’s pairwise comparison matrix (see Table 3). Euclidean distance 
process and reclassification of images was carried out and suitability images finally derived using 
Equation | based on AHP and Binary methods. The EDFs were finally used as independent variables 
in a multiple regression analysis as shown in Equations 2 and 3 for change otherwise known as 
sensitivity analysis. This process allows the analyst to identify the less important factors, thereby 
making use of the most relevant factors. 


Table 2: Lists of Suspected External Driving Factors (EDFs) 


Geophysical Slope SRTM 1 arcsecond 30m GSD 
Elevation 
Soil 


Neighbourhood effect Land use Land cover Pixel 
Proximity to environmentally | Distance to market/Commercial area | Euclidean distance to nearest 
sensitive/prominent areas Distance to industrial area (Km) 
Distance to water body 
Distance to mines area 
Distance to river 
Channels 
Distance to settlements 
Distance to conflict areas 
Proximity to Utilities Distance to Major roads Pixel 
Distance to railway Euclidean distance to nearest 
Distance to water facility (Km) 
Distance to Health care 
Distance to Educational area 





Table 3: Adopted thresholds for suspected EDFs 


Factors Distance 
Major road 
Water facilities 
Health care 
Education 
Market/commercial 
Industrial 


SWAN NB WN Re 


ioe) 


Water body 100.00m 
Mining area 100.00m 
9 Railway 50.00m 


Slope <5" 
River 30.00m 
Conflict area 1000.00m 
Settlements 1000.00m 





2.3.2. Generate suitability images 

Euclidean distance images were first generated for all the EDFs followed by reclassification process 
based on the thresholds adopted in this work as shown in Table 3. The reclassified images based on 
AHP and Binary methods were multiplied through raster calculator in ArcGIS 10.1 by the 
corresponding weights computed as shown in Table 4. Therefore a total of thirty two (32) images 
were generated with sixteen each from AHP and Binary methods. Equation | shows how each of the 
suitability index maps was created. 


Si; = Xij-Wn- Cn (1) 


Zitta et al., 2022 61 


Nigerian Journal of Environmental Sciences and Technology (NIJEST) Vol 6, No. 1 March 2022, pp 58 - 70 


Where: 
Si = Suitability indexes for cell ij at time t 
Xj, = Criteria m at cell ij at time t 


Cy, = Constraint value 
W,, = Weight for criteriam 


Suitability index map Si : base on EDFs generated through AHP model is given as U; j_AHP while 
EDFs generated using Binary model is given as Vij p m- it is important to state that all the factors do 
not influence the rate of urban growth in the same way due to polarity. For instance, a high population 
density of a region will create more pressure for the conversion of natural resources. Contrary to this, 
the rate of urban growth will be high in those areas near the city centre. Therefore, a fuzzy set 
standardization procedure was followed to remove the conflicting character in the dataset and to make 
all factors unidirectional with an equal scale. Images were standardized to a common scale of 0-255 
that is all the factors were brought to a common measurement unit or scale (suitability). A value of 
255 indicates the highest suitability and a value 0 indicates the lowest suitability of that particular 
category. 


2.3.3. Multiple regression analysis (MRA) 

The general multiple regression (MR) formula as shown in Equation 2 was performed to analyze the 
relationship of one or more independent variables on a single dependent variable. The regression was 
conducted for image files of the various EDFs (independent variables) and the LULCC of t, — t3 as 
the dependent variable. This was performed in Idrisi selva 17.0 software through GIS 
analysis/statistics/multireg module for both AHP and Binary suitability images. The Multiple 
regression equation for AHP and Binary suitability images takes the form as is given in Equation 3 
and 4 respectively. 


Y=A+Bx,+ Cx2 + Dx3 (2) 
Yat b,U;j ane, + b2U;j np, aie by Uj j_anP, (3) 
Where: 


Y = dependent variable (LULCC from t, — f3) 

a = constant 

b, — Coefficient 

Vij p u,& U; j_AHP 4 _ Independent variable (suspected external driving factors for the 
period t; — t3 1.e t= 33 years) 

n = Number of independent variables (16) 

t;= LULC of 1986 

t,= LULC of 2000 

t3= LULC of 2019 


2.3.4. Final selection of appropriate external driving factors (EDFs) 

This was achieved through the result of MR by considering the corresponding sign either positive or 
negative as displayed by regression coefficients (r) which is also known as slope coefficient. 
Independent variables with consistent behavior from both models were selected as the final 
appropriate EDFs. A positive sign means as the dependent variable increases, the independent 
variable also increases and both variables move in the same direction while a negative sign indicates 
an increase in the dependent brings about a decrease in the independent variable and both variables 
move in opposite direction (Jim, 2019) . 


3.0. Results and Discussion 


The results showing weights of all the EDFs are displayed in Table 4 after formulation of pairwise 
comparison matrix. A consistency ratio (CR) of 0.017 was achieved which is slightly above the 0.1 
target, and the weight sum up to 1.0000. It can also be deduced that, factors like Elevation, distance to 
market and distance to river are having higher weights while slope, soil, distance to mining area and 
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distance to water body are having the lowest weights. This may be attributed to one of the major 
drawbacks of method adopted (Saaty’s pairwise) that is, varying expert’s opinion. Table 5 is the result 
of MR performed on the thirty two (32) EDFs with sixteen (16) each from both methods after 
performing fuzzy set standardization. The analysis of this result shows that LULC, distance to conflict 
areas, distance to river/channels and distance to settlements are proximity to environmentally 
sensitive/prominent areas while distance to major roads is proximity to utilities. 


Table 4: Computed weights based on Saaty’s pairwise comparison matrix 


EDF 


CI = 0.27164 RI = 1.5978 n= 16 Amax = 20.07459 
CR = 0.17000 





Table 5: Selected appropriate EDFs 


r coefficient (Binary) r coefficient (AHP) 
0.001417 0.183332 
LULC 0.016682 2.702839 


0.001730 0.040826 
0.005624 0.051394 
Distance to river| channels | 0.003778 0.000784 





Result in Table 5 above shows the EDFs that were consistent in their influence or behaviour against 
the change period of study. In this study, emphasis was on the sign (direction of the relationship) not 
the magnitude of coefficient as the case has always been. Direction of relationship is very important in 
this research because, much attention is drawn to spatial movements. As landscape changes to which 
direction, what are the influencing EDFs that are responsible for those changes in that direction? The 
images as shown in Figures 2a — 2e are the results derived from AHP method, while Figures 2f — 2] 
are from Binary method. These are the final selected EDFs contributing to the changing landscape of 
the study area. These images shows areas that are suitable and not suitable for urban land 
development in the study area as presented using two different methods. The AHP image displays a 
suitability range from 0 — 255, that is suitability of any portion in the study area is dependent on the 
pixel value. Similarly, for binary images, they are | and 0. Areas with pixel value 1 are suitable while 
areas with 0 are not suitable for urban land development. 
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Figure 2a: Distance to conflict areas Figure 2b: LULC 
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Figure 2c: Distance to major roads Figure 2d: Distance to river channel 
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Figure 2e: Distance to settlements Figure 2f: Distance to conflict areas 
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Figure 2g: LULC Figure 2h: Distance to major roads 
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Figure 21: Distance to river channels Figure 2j: Distance to settlements 


Distance to major roads: the influence on landscape change is predominantly in the north running 
through the central towards the east and down to the southern part of the research area as seen in 
Figure 2a & 2h. In Adzandeh et al. (2015), the rate of land development in Jos is alarming. This is 
further confirmed in this research as a result of massive roads construction, other sub-urban 
settlements like Kwang, Lamingo Haske, Gwafang, Ten Commandments, Guratop and Du have 
witnessed increased urban land development. These newly constructed roads have greatly influenced 
the cost of landed properties and standard of living in these areas. They are now major alternative 
routes (bye-pass) for transporters easing traffic gridlock usually experienced in the city centre. Land 
use Land cover (LULC): as seen in Figures 2b & 2g respectively, there is little or no effect of 
potential land development in the northern part through the central down to the south of the study 
area. This is due to the existing developments, but the LULC is a pointer to other parts like the west 
and some part of south east as highly suitable for developments. 


Distance to river channels: apart from the main river channels as shown in Figures 2d & 21, other 
areas around the EDF clearly indicates the suitability of land development in all directions of the 
study area. This is quite typical of the settlement expansion witnessed today in Jos as most residents 
along this EDF takes advantage of proximity and engage in both dry and rainy season farming. 
Distance to conflict areas is also another influencing factor of land development as shown in Figures 
2a & 2f respectively. The dark spots on the images are vulnerable areas and can be seen around the 
south, some part of east and south west areas. Suitability of land development in the study area is 
towards the central to the northern area and some parts south east. Also, Adzandeh (2015) pointed out 
that these areas are constantly witnessing a decline in urban land development since 1984. Distance to 
settlements: the spatial distribution of existing settlement has equally influenced the development or 
existence of other settlements due to proximity. Figures 2e & 2; shows that almost all parts of the 
study area are suitable for urban land development as a result of the already existing settlements. Even 
though, in the Binary suitability image it shows that only a small portion in the south east is not 
suitable. This implies that, it is safer to live close to each other than stay in isolation because of 
prevailing factors like conflicts and insecurity. 


4.0. Conclusions 


To determine the drivers of change, sixteen (16) suspected external driving factors (EDFs) were 
identified as those of geophysical, neighbourhood effect as proximity to environmentally sensitive 
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areas and proximity to utilities. Euclidean distance analysis was carried out on the proximity EDFs 
and reclassification based on AHP scale | — 9 and binary mode as | & O. The Saaty’s pairwise 
comparison matrix was generated and weights for all the EDFs were computed with a consistency 
index (CI) of 0.17. Suitability images were generated by multiplying each EDF with the 
corresponding weight. The fuzzy set membership standardization was carried out on the AHP image 
through sigmoidal function type with monotonically increasing membership function. The Multiple 
regression technique was used to measure the sensitivity of all the EDFs (images) against the change 
image (1986 - 2019). Based on the regression analysis, five (5) EDFs were consistent in both 
approaches (AHP and Binary). They are: LULC, distance to conflict areas, distance to river/channels 
and distance to settlements are proximity to environmentally sensitive/prominent areas while distance 
to major roads is proximity to utilities. 
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